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ABSTRACT 

The statistics of giant arcs and large separation lensed quasars provide pow- 
erful constraints for the parameters of the underlying cosmological model. So far, 
most investigations have been carried out using pure dark matter simulations. 
Here we present a recipe for including the effects of baryon cooling (i.e. large 
galaxy formation) in dark matter A^-body simulations that is consistent with ob- 
servations of massive galaxies. Then we quantitatively compare lensing with and 
without applying this baryon correction to the pure dark matter case. Including 
the baryon correction significantly increases the frequency of giant arcs and lensed 
quasars, particularly on scales of 10 arcsec and smaller: the overall frequency of 
multiple images increases by about 25% for source redshifts between Zg = 1.5 
and 7.5 and splittings larger than about 3 arcsec. The baryon rearrangement 
also slightly increases the fraction of quadruple images over doubles. 

Subject headings: cosmology: gravitational lensing, arcs, quasars, galaxy clusters 



Introduction 



It has been realized for some time that the abund ance of clusters of galaxies can provide 



very stringent constraints on cosmological models (e.g. Henry fc Arnaudlll991 : Oukbir fc Blanchard 
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1992 



1996 



White. Efstathiou. fc Frenklll993l : iBahcall fc Cenlll993l : lEke. Cole fc Frenklll996l : IViana fc Liddle 
Bahcall fc I2OO3I . and references therein). The reason for this is that — since only 



5% to 10% of the stellar mass is in rich clusters, the exact fraction being fixed by the def- 
inition of "rich" — they represent relatively rare fluctuations of approximately two sigma. 
Consequently, as they lie on the tail of the density distribution, relatively small changes 
in the mean amplitude of the cosmic perturbation spectrum can produce large changes in 
the expected numbers of rich clusters. In fact, gravitational lensing by clusters, which de- 



pends on the most dense fra ction o 



of cosmological parameters ( 


Turner 


1990 


Li & Ostriker 


2002. 


2003: 


Li et al. 


2007 


)• 



these systems, provides an especially sensitive measure 



in the normalization of the fluctuation power spectrum. For example, Spergel et al. (2006) 
give an estimate based on WMAP 3-year and SDSS data that is significantly lower than the 
WMAP 1-year data, but Evrard et al. (2007) argue that the cluster X-ray data indicate a 
higher normalization cosmology. It is hoped that efforts such as the present paper, based on 
ray-tracing through specific cosmological models, can help resolve this controversy. 

The frequency of giant arcs and large separation quasar lenses has received increasing 
interest in the last few years as a powerful tool for constraining cosmological parameters. 
The recent discovery of a few widely separated qu a sar le i ises (A^ > 10 arcsec) from the 



Sloan Digital Sky Survey fllnada et al.l l2003l . l2005l . I2OO6I : lOguri et al.l l2004l : ISharon et al. 



20051 ) has also spurred such studies. Almost all theoretical studies using iV-body simu 



lations are constrained on dark matter onlv (Barte 


mann et al. 1998: 


Wambseanss et al. 


1995 


19981: Meneehetti et al.l 20031: 


Wam 


3se:anss et al. 


2004: 


Dalai. I 


lolder & Hennawi 20041: 


Oeuri & Keeton 20041: 


Ho & White 


2OO5I; 


Horesh et al. 


2OO5I: 


Li et al. 


2OO5I 


20071: Hennawi et al. 


2007a 


b: Hilbert et al. 


2007). In such studies it is usuallv assumed that the barvonic matter 



essentially follows the dark matter particles. This may be true in general on large scales, 
but it certainly is not a very good assumption in the central parts of galaxy clusters or for 
very massive isolated galaxies, which are precisely the essential matter concentrations for 
the production of giant arcs and wide quasar lenses. The baryonic component is also known 
to dominate in the central regions of elliptical galaxies. 

Recent studies have begun investigating the influence of baryons on DM clustering. 
Puchwein et al. (2005) looked into the effect of gas physics on strong lensing by individual 
galaxy clusters. They found that cooling and star formation can increase the strong-lensing 
efficiency considerably. Jing et al. (2006) studied the influence of baryons on the clustering of 
matter and weak-lensing surveys, and found that the clustering of total matter is suppressed 
by about 1% on large scales of l/iMpc~^ ^ k < 10/iMpc~^, while it is boosted between 2% 
and 10% on small scales of ~ 20/iMpc~^; they conclude that this should be measurable 
with future weak lensing surveys (Jing et al. 2006). Lin et al. (2006) looked into the 
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influence of baryons on tlie mass distribution of dark matter halos and found an increase of 
the concentration parameters by 3% to 10% as compared to pure A^-body simulations. Rozo 
et al. (2006) studied the effect of baryonic cooling on giant arc abundances for individual 
clusters; they found that the arc abundances can be increased by factors of a few. 

In this study we quantitatively investigate the effect of baryon cooling on the statis- 
tics of arcs and widely separated multiple quasars. We first motivate and describe a simple 
recipe for redistributing part of the matter in the centers of halos to approximate the effects 
of cooling and star formation, and present some tests of this method. Subsequently we sum- 
marize our use of the ray shooting method. Then we describe the quantitative results of our 
baryon redistribution recipe by directly comparing pseudo-three-dimensional ray shooting 
simulations with and without this baryon cooling, and finally discuss this effect with respect 
to the observational situation. 



2. Methods 

2.1. Baryon rearrangement 

In order to approximate the effects of galaxy formation in our ray shooting simulations, 
we locate halos in all lens planes and identify the amount of mass likely to have cooled into 



Such an approach is supported by observations, e.g. 



Pens; et al. 


( 


2004 


) or 



( 120071 ). who have found that early- type galaxies are consistent with an isothermal profile 
out to > 300 h~^kpc from both kinematical and lensing data. Here we describe in detail 
the ways this rearrangement is performed. The approach that we adopt is admittedly very 
rough, and so, as in all semi-analytical modeling, we constrain the free parameters to respect 
observational data. The requirements the model must fulfill are: a) The fraction of the 
baryonic component which is rearranged into stellar systems is consistent with the stellar 
mass fraction of the universe (cf. Figures [2] and [3]) . b) The cosmic buildup of mass in stellar 
systems parallels what is known from observations (Figure [3]) . c) Individual mass profiles 
within the baryonically dominated systems are consistent with kinematic and gravitational 
lensing data (Figure H]). d) The distribution of systems as a function of stellar mass and 
circular velocity approximates observational data (Figures [2] and H]). 

Suppose we have projected the mass in some volume onto a plane for the purpose of ray 
tracing (as described in the next section). Let us denote the surface density of a given pixel 
as a and the pixel size as /, so that the mass in each pixel is cx/^. For a given 2-D plane, each 
pixel is examined in turn to see if it has a higher surface density than any other pixel within 
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a radius roughly the size of the smallest objects resolved in the simulation, Rb = 30h~\pc 
physical. If so, this pixel is taken to be the center of a halo. For a given cooling radius Rc (to 
be determined below), the background surface density a^g is set to the mean surface density 
between R^ and Rc + {Rb/2). The mass of the halo is defined as 

Mh = j W{R/Rc) [a(R) - a^g] dA , (1) 

integrating from the halo center out to Rc, here we will use the window function W{x) = 
[1 — x^]^. With the mass weighted radius of the halo given by 

Rf, = j RW{R/Rc) [a{K) - a^g] dA , (2) 

the halo temperature is defined as 

Gmp\ Mh 



with G and k being the Newton and Boltzmann constants, respectively, and mp the proton 
mass. 

We wish to find a cooling radius such that half the baryons contained within this ra- 
dius will have formed stars. For ionized gas in a galaxy clust er the coo ling time from 



Bremsstrahlung radiation is given by t^. = 9 x 10^Tg^'^^?T,~^yr (cf. ICoxll2000l ): here rif. is the 
electron density per cm^, which we will take to be ng = ?)Mh/ {4:TiR\mp). Let Tq be the 
temperature at which the cooling time equals the age of the universe at the redshift under 
consideration, thuh'- 

2 / ^hub \ ^ „87 



To set Rci we start with a small estimate for the cooling radius, and increase it until 

e-^-'^^^ = h , (5) 

with the cooled fraction = i, is satisfied. A number of oversimplifications have gone into 
the calculation of and Tq, so the variable t] is introduced to account for these approxima- 
tions. This parameter can be adjusted to give more or less stellar mass; we find = \ gives a 
reasonable amount of stellar mass, consistent with observational constraints. Problems will 
arise with this algorithm if the aperture Rc increases to the point where it begins to include 
neighboring structures, as can happen in the case of a galaxy-sized halo inside a group or 
cluster. Thus, if during this process it happens that Rh/Rc > ^, we stop increasing Rc and 
leave larger than a half. 
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Mass representing the cooled baryonic component is then removed from inside i?c, re- 
ducing the surface density to a temporary value 

a,(R) = a{R) - ^f,W{R/R,) [(t(R) - a,g] , (6) 

with the caveat that any pixels denser than the central pixel are left unchanged. This removed 
mass is then added back in, using an isothermal profile A/R within a smaller radius Ra- 
With the isothermal component added back in, the final surface density is 

af{R) = a,g + AW{R/RA)/R • (7) 

For a given Ra (set below) the value of the constant A is set by conservation of mass, i.e. 
the mass put back in, Mg, must of course equal that taken out originally. Figure [T] illustrates 
how our method of baryon redistribution affects the velocity profile of a halo. The four lines 
in this example halo show the initial circular velocity from the pure N-body simulation (the 
solid line), the circular velocity after mass removal as given by Eq. [6] (dotted), the velocity 
curve of the isothermal component to be added back in according to Eq. [7] (dot-dashed), 
and the final circular velocity profile after the redistribution is complete (dashed line). This 
example corresponds to a brightest cluster galaxy residing at the center of the cluster dark 
matter halo. 

Figured shows the "stellar mass function" (that is, the distribution of Mg) in halos 
identified by our algorithm at the redshifts z = 0.05, 0, 55, 1.05. The distribution of stellar 
masses is quite similar at all three redshif ts. The lower dotte d line results from taking 



the luminosity fun ction of SPSS galax ies (IBlanton et al.l l2003l ) and assuming a constant 



Mgtar /L = A A (see iGayazzi et al.l 120071 ): the upper dotted line instead assumes Mstar/L oc 
L'''^^ ( IVale fc Ostrikerl I2OO6I ). As a further test, summing the total mass rearranged at a 
given redshift yields an estimate of the stellar mass in large spheroidal galaxies. The points 
in Figure [3] show the density of this component (as a fraction of the critical density) versus 
redshift. The error bars give the error on the mean value of the 486 planes used; this 
error becomes larger at low redshift primarily because the size the planes is smaller. The 
total fraction of mass that is rearranged rises slowly from VLgtar = 0.7 x 10~^ at 2 = 1 to 
^star = 1-0 X 10"'^ aX z = 0.29, remaini ng at this value for lo wer redshifts. Also shown are 



empirically based estimates (compiled in lNagamine et al.ll2006l ) for the total mass in elliptical 



galaxies. The solid line shows the total mass in elliptical ("spheroidal") components, and 
the dotted line indicates an estimate of the mass limited to systems with > 215km/s 
which would produce a splitting angle of >3 arcseconds. These hues bracket our points, 
indicating that we are not overproducing galaxies, but are including those hkely to produce 
significant lensing. From these tests we conclude that our choice of ?], which determines the 
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mass rearranged in a given halo, is reasonable. The stellar mass function is insensitive to this 
parameter — doubling rj causes very little change to the mass function above 2 x lO^^h'^MQ. 

The v a lue o f Ra is set to match the o bserved size distr i bution of SDSS galaxies. 
Shen et al.l (120031 ). using the stellar masses of iKauffmann et al.l (120031 ). find that the Ser- 
sic half-light radius of early type galaxies varies with stellar mass as M^'^^. When using the 
above procedure on the lowest redshift planes, we find that the cooling radius Rc varies with 
the mass removed as Rc oc M^'^°. Thus by setting 



Ra — Rc 



2.4 X W^h-^Mr. 



0.36 



the half-mass radius of t he added mass ma tches very well the half-light versus radius relation 
for early-type galaxies of lShen et al.l (120031 ): this is shown at low redshifts in the bottom panel 
of Figure |H If this choice of Ra would result in mass being added to a pixel which is below 
the background density, i.e. cxj < abg, we instead reduce Ra to the level where this will not 
occur; it is quite rare that this adjustment is invoked. 

The value of Ra also sets the isothermal circular velocity Vc = 2^/GA. The upper panel 
of Figure m displays the circular velocity measured at Ra as a function of stellar mass. These 
velocities agree reasonably well wi th the relatiori between stellar mass and rotation velocity 
at 2.2 disk lengths measured by jPizagno et al.l (120051 ) . who find that V2.2 increases with 

1 /3 

stellar mass roughly as Ms . At higher stellar masses our circular velocity is about 10% 
higher than the measured ¥2,2- A lower limit on the isothermal circular velocity is imposed: 
for Vc < 125km/s, the halo is simply returned to its original density profile. This velocity 
cut also puts a lower hmit on the splitting angles affected by this procedure, which is of 
order one arcsecond. 

To summarize this section, the simple procedure described here reproduces the stellar 
mass function of large galaxies as a function of redshift, with the appropriate sizes and 
circular velocities. In Section 3 we investigate what effects these galaxies with rearranged 
matter distribution have quantitatively on the lensing frequency. 



2.2. Mass planes, ray shooting and identification of arcs 



The methods used to create the mass dis tribution in the len s ing p lanes and to carry 
out the ray tracing are described in detail in IWambsganss et al.l (120041 ). Here we briefly 
summarize the procedure and give the numerical values of our parameters. Essentially, light 
rays are followed backward from the observer through a series of lens planes, approximating a 
three-dimensional matter distribution; the lens planes are obtained from A^-body simulations. 
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The parameters of the underlying cosmological model are: Qm — 0.3, Qa — 0.7, Hq — 70 

km/sec/Mpc, (78 =0.95, and ns=l. The iV-body simulation was performed in a box with a 
comoving side length of L = 320/i~^Mpc. We used N = 1024^ particles, so the individual 
particle mass is rrip = 2.54 x lO^h^^ Mq. The cubic spline softening length for all particles 
was set to e = 3.2h^^ kpc. The output was stored at 19 redshift values out to z ^ 6.4, such 
that the centers of the saved boxes matched comoving distances of (160 + nx 320)/i-^Mpc, 
where n = 0, 18. 

For 243 lines-of-sight traversing the box, two lens planes were produced by bisecting 
along the line-of-sight and projecting the mass in each 160/i~^Mpc-long volume onto a plane. 
In the lowest redshift box the Icnsing planes are 1.9/i^^Mpc on a side, and the pixel size is 
I — 2.3/i~^kpc. With increasing lens redshift, we keep the number of pixels constant but 
increase the physical size of the planes, such that the pixel opening angle remains constant. 

Light rays are then propagated backwards through these lens planes, beginning with a 
regular grid at the lowest redshift lens plane (i.e. the image plane) and working to higher 
redshift by considering the proper deflection in each lens plane. For each source redshift, 
the coordinates of the rays in the image plane and in the source plane are stored. For 
analyzing the imaging properties, we impose a regular grid of sources in the source plane 
and identify for each source position the image multiplicity (by far most of them are single 
images, some are triples, we also identify quintuples, septuples etc.), and then for each image 
the corresponding image position and image magnification. 

In the hterature, different authors use different definitions for what an "arc" is. As stated 
in Wambsganss et al. (2004), we demand as a miminum requirement that an arc be part of a 
multiple image system and have a certain minimum magnification (we chose different values 
between 5 and 25 as threshholds). Since strongly lensed sources are almost always highly 
distorted in the tangential direction, in particular the strongly magnified multiply imaged 
ones, the magnification is in general close to the length-to-width ratio of such an image. We 
double-checked this assumption in a number of individual cases and found good agreement 
above the 90% level. 

For the multiply imaged sources we order the images by magnification and determine 
the distances between the images, whereby the "image separation" is defined as the angular 
distance between the brightest and the third brightest image. More details are given in 
Wambsganss et al. (2004, 2005). 
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3. Ray shooting Results and comparison with observations 

3.1. Ray shooting 

We performed two complete sets of rayshooting simulations as detailed in Subsection 
2.2, once with the original pure A^-body matter screens and once with the new recipe of 
rearranged matter, creating 100 independent realizations for each of the two. This procedure 
allowed us to study the effect both statistically and on individual lines-of-sight (the latter 
turned out to be essential for a complete understanding of the results, as described below). 
Three source redshifts are considered, Zg = 1.5, 3.7, and 7.5. 

With the rearrangement, the number of multiple images increased for all source redshifts, 
particularly for angular separations up to 15 arcsec. This is shown in Figure [5l where the 
multiple image cases are binned according to the image splitting using bin widths of 5 arcsec 
(note we are highly incomplete in the first bin, because we do not properly resolve image 
splittings smaller than about 3 arcseconds). The level of increase can be read off of the 
integrated distributions in Figure [6], where the dashed lines show the integrated frequency 
of multiple images as a function of separation for the pure N-body case, and the solid lines 
show the scenario with baryon rearrangement. Qualitatively, it is obvious that including the 
baryon redistribution produces more multiple images. The amount by which the rearranged 
case is higher than the pure A^-body case is almost independent of source redshift: the 
fraction is about 25% for the redshift range Zg = 1.5 to 7.5 (it varies between 22% and 28%). 
Beyond the total lensing frequencies, it is interesting to study the differential e^iect: it is for 
relatively small separations (< 10 arcmin) in particular that the baryon-redistribution case 
leads to more multiple images: about a 70% increase in the lowest bin, and 30% increase 
in the 5-10 arcsecond bin. This excess decreases with increasing separation, such that at 
a splitting of ~ 12 to 15 arcsec the two cases result in the same frequencies for multiple 
images. This effect is easy to understand: the baryon redistribution steepens the innermost 
parts of the density profiles of halos, and hence can bring DM halos which are originally just 
below the critical density for lensing above this value, thus producing multiple images. 

At still larger separations, the baryon redistribution case apparently produces slightly 
fewer cases (at the few percent level). The reason for this is not quite so obvious; in fact, 
it is rather counter-intuitive: if a fraction of the total mass becomes more concentrated, the 
total mass inside an image pair should be the same or larger. The explanation is that this 
effect is just an artifact of our definition of image separation. As described above, we define 
the separation of a multiple image system as the angular distance between the brightest 
and the third-brightest image. Our method always identifies the "odd" image in triples and 
quintuples, although in practice most would be identified as double and quadruple images. 
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because the odd images in most observational situations are demagnified and not detectable 
(hence we refer subsequently to those CclSGS clS "doubles" and " quadruples" ) . A detailed look 
at the distribution of multiple images shows that a side effect of the redistribution of baryons 
is that the number of quadruple images increases more strongly than the number of doubles 
(see Figure [7]). We are able to verify this by comparing individual lines of sight, with and 
without the matter rearrangement. When we compare one-to-one a situation that produces 
a double image in the DM-only scenario and a quadruple image in the baryon-redistribution 
case, the additional image pair occasionally is brighter than the faint "odd" image, and it 
happens that this pair is close to the brightest image of the multiplet, resulting in a smaller 
measured separation than in the pure N-body case. This leads to the apparent small deficit 
of large separation cases in the redistributed-baryon case. 



3.2. Comparison with observations 



Over the last decade, v arious papers h a ve inv e stigated the obs e rvatio nal oc currence of 



giant arcs statistically, e.g. iLuppino et al.l (Il999l ). iGladders et al.l (120031 ) . and ISand et al. 
( 120051 ). The frequency of arcs reported in these studies varies from one giant arc per 45 
deg^ to about one arc per 10 deg^; the large spread may be explained partly by the slightly 
different definition and limiting magnitude, and may be related to the relatively small number 
of cases per study. Overall, just about 50 arc systems had been used for these statistical 
analyses, found in various surveys with diverse selection criteria. 

Two recent studies b ased on SDSS data present quite different results as well. On 
one hand, lEstrada et al.l (120071 ) report no arcs found in a systematic inves tigation of 825 



SPSS clusters (other than one serendipitous discovery). On the other hand, iHennawi et al. 



(120061 ) presented first results from a new survey for giant arcs, which may as much as double 
the number of known arcs (about 30 new systems already reported). This study is very 
promising not just due to the large number of new arc systems to be expected, but - at least 
as importantly - also by the fact that it will find them by well defined selection criteria. 

Until these future results will be published, we c a n com pare our results here (although 
only "differentially" ) to those in IWambsganss et al.l (120041 ) . There it was found that the 
simulated pure N-body based arc frequencies are consistent with the observations, if we 
allow the source redshifts to extend well beyond z > 1 {as the distribution of observed arcs 
is as well). The result presented here — that by including the baryon rearrangement the 
predicted number of arcs increases by about 25% — puts the predicted arc frequency a bit 
on the high side, though still with considerable uncertainty both on the observational (e.g., 
absolute numbers, selection criteria) and the theoretical side (e.g., parameters of underlying 
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cosmological model, exact definition of "arc"). Furtliermore, tlie model used liere lias a 
high normalization for the amplitude of fluctuations as compared to the recent res ults from 



WMA P; a lower normalization would significantly reduce the lensing frequency (ILi et al. 
2007h . 



4. Summary and Conclusion 

This paper quantitatively investigates the effects on lensing of baryon redistribution in 
originally pure dark matter A^-body simulations. We describe a heuristic recipe for rearrang- 
ing the baryons in dense dark matter halos that is consistent with observational data for 
massive galaxies, present a number of tests of this prescription, and then apply it to a cos- 
mological simulation used for multiple lens plane ray shooting. We compare the frequency 
of multiple images, the image separation and the image multiplicity between the original 
A/"-body matter distribution and the redistributed version. We find that on average the case 
taking into account the redistribution of baryons produces 25% more multiple images; this is 
almost independent of the source redshift in the range Zs = 1.5 to 7.5. For splittings between 
5-10 arcseconds, the number of multiple images increases by roughly 30%, and by 70% for 
smaller separations < 5 arcsecond. But this last result must be treated with caution since 
we are resolution limited for splittings < 3 arcsecond. As noted, most of the new multiple 
images systems occur for angular scales < 10 arcseconds. We also find that the number 
of quadruple images increases more than the number of double images. This produces an 
apparent slight reduction of larger separation cases, which is in fact only an artifact of the 
way we define image separation, namely as the distance between the brightest and the third- 
brightest image. Since arc statistics are such a good tool for distinguishing between models 
with different cosmological parameters, we would like to emphasize once again the need for 
both very good observational studies and for more explorations of the important parameters 
in the simulations. 
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Fig. 1. — The circular velocity profile of an example halo, demonstrating the effect of redis- 
tribution. Solid line: initial circular velocity. Dotted line: after mass removal (cf. Eqn. [6]). 
Dot-dashed line: velocity of the isothermal component (cf. Eqn. [7]). Dashed line: the final 
circular velocity profile. 
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Fig. 2. — Solid lines: stellar mass function as determined with the recipe described in the 
text, at redshifts ^=0.05, 0.55, and 1.05 (top to bottom). The lower dotted line is taken 
from the SDSS luminosity function (Blanton et al. 2003, Gavazzi et al. 2007), assuming 
Mstar/L — 4:A; the upper dotted line instead assumes Mstar/L oc L°-^^. 
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Fig. 3. — The total rearranged "stellar" mass as a function of redshift (determined by 
applying the baryon redistribution recipe described in Section 2.1 ). The solid lirie is a n 
empirically based estimate of the total mass in elliptical galaxies from lNagamine et al.l (120061 ): 
the dashed line is limited to systems with Vc > 215km/s, which would produce a splitting 
angle of >3" . 
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Fig. 4. — Top panel: Circular velocity V at Ra as a function of stellar mass M^teiiar-, for 
halos identified in the ^=0.05 and 0.17 lens p lanes. The line is the circular velocity measured 
at 2.2 disk lengths by iPizagno et al. Bottom panel: Half- mass radius R1/2 of the 

redistributed mass as a function of stellar m ass Mf-tpUn r - Th e line is the Sersic half-light 



radius of early type galaxies from SDSS, see IShen et al.l (120031 ) 
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Fig. 5. — Frequency of multiple images as a function of separation (binned with 5 arcsec 
width). The number of multiple images clearly increases with decreasing splitting angle; 
the drop in the smallest separation bin is due to the fact that we are highly incomplete for 
splitting angles below 3 arcseconds. The N-body case is shown in the top panel, the sce- 
nario with baryon redistribution in the bottom panel. The histograms reflect three different 
source redshifts: = 7.5 (highest, black, loose shading), = 3.7 (middle, red, intermediate 
shading), z^ = 1.5 (lowest, green, densest shading). In the baryon redistribution scenario 
we find more multiple images (on average about 25%); this is most pronounced in the small 
separation regime. 
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Fig. 6. — Integrated probability distribution of multiple images as a function of separation 
for a dark matter only simulation (green dashed line) and with baryon redistribution as 
described in the text (blue solid line). Top panel: For a source redshift of = 1.5; middle 
panel: = 3.7; bottom panel: Zs = 7.5. 
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Fig. 7. — Image multiplicity as a function of source redshift for "triple" images (triangles) 
and "quintuple" images (pentagons). The pure N-body case is shown with solid symbols 
and dashed lines, for the baryon-rearrangement case we used open symbols and solid lines. 



